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SIWOPSIS 


A non linear mathematical model for the turbine 
of the 110 MW unit of Panld. Thermal Power Station is deve- 
loped. The model for the unit consisting of the high 
pressure, medium pressure and low pressure turbines, the 
reheaters, the condenser and the high pressure and low 
pressure feed water heaters, is simulated on a digital 
computer and open loop responses for a 10 MV/ load drop and 
step change in heat input to the rehcaters have been 
obtained. Various results arc discussed. Plant test 
conducted for a 10 MW load drop is compared with the model 
but extensive comparison of model responses with plant 
tests could not be made because of the lack of proper 
instrumentation available in the unit. A review of litera- 
ture has been done and recommendations for further work are 
given . 
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CHAPTER 1 


INTRODUCTION 

Economic prosperity and social progress of a 
developing countr3^ like ours is heavily dependent on how 
we utilize our power generation resources. There are 
three major sources of energy viz. Hydrc_ power, Thermal 
power and Nuclear power. The hydro power potential of our 
country is large hut it is restricted hy factors like large 
initial investment, long gestation period, excessive depen- 
dence on monsoon etc. Nuclear power generation prospects 
are not very bright due to economic, political and tech- 
nological considerations. As a result, the current trend is 
to go for more and more fossil fuel fired plants. Of the 
existing generation capacity of 26,000 MW, the hydel sources 
account for 4-0 percent. The next five years plan projec- 
tion envj. sages augmentation of Installed capacity by 
18,500 MW. Of this 13,000 MW will come from coal fired 
thermal stations, 3,800 is assigned to hydel and the 
balance to Nuclear plants. So bulk of the power will come 
from themal and super thermal power stations. 

A modern power station is a gigantic complex 
by itself and consists of several sub-systems and innume- 
rable components. The performance of the overall plant 
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depends critically upon the functioning of these components 
and sub- systems. Besides, demand for higher efficiency is 
leading towards higher steam pressures and temperatures. 
There is also a tendency to go in for larger and larger 
unit sizes for reduction of capital costs. On the other 
hand, more and more stringent requirements are imposed 
regarding safety and control. All these have 3 -ed to the 
development of modern power systems and power stations 
with highly complex instrumentation and control systems ^ 
quite often aided 137- computers for on line data logging, 
monitoring and in a few cases for control purposes. 

It is in this context that understanding of the 
power plant dynamics becomes important. Smooth function- 
ing of the plant depends upon the nature of the interaction 
between the various subsystems. How well a plant behaves, 
is determined not only by its steady state performance, 
but to a large extent by how it reacts to disturbances and 
abnormal conditions. These disturbance consist of momen- 
tary or long time malfunctioning of internal components or 
external disturbances. Plant control and safety system 
must override such situations and outages. Effective 
steps can only be taken if behaviour of the system is 
fully understood and only then can suitable control and 
safety systems be engineered to meet the situation. 
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Dynamic simulation in such conditions becomes a 
very important tool because of folloi-n.ng reasons ; 

1. It is not feasible to create major power system distur- 
bances for the purpose of examiiiing system behaviour 
and effectiveness of proposed strategies or actions. 
Experimental evaluation of strategies is therefore 
impractical. Simulation can be used to judge the 
merits of these strategics under such major disturbances. 
Moreover sensitivity of these strategies to system 
changes and system model parameters can also be deter- 
mined. 

2. By developing models for different components much 
insight is gained into the unit dynamic characteristics, 
control potential and system behaviour, 

3. Simulation provides a means of evaluating the system 
behaviour and helps in predetermining the critical 
parameters and their response to operator action etc. 

Power system dynamics (upto several minutes following 
a major system upset) are not well understood and 
simulation helps in this regard. 

Generally three lands of problems are studied 
in dynamic analysis of power systems. First, the study 
of instability that may occur when there is a fault on 
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the generator side. This kind of study is done over 
few seconds of the plant operation. Second, the study 
of the response of the system to load changes. This study 
is done over few minutes. Third, the study of start 
up and shut down etc. which are carried over for few hours. 

It is the second category in which the present 
problem falls. The first step towards conducting such a 
study is developing a mathematical model which adequately 
describes the system behavioLir over the interval needed for 
that particular problem. This requires a thorough under- 
standing of tho component systems, and processes and their 
interaction. 

As the features of a particular unit differ 
from one power plant to another, need was felt for develop- 
ing a model for a unit which is most representative of the 
current boiler turbine units being used in our country. 
Proximity to I.I.T. Kanpur weighed heavily in selecting the 
110 MVJ unit of Panki Thermal Power Station for this puipose. 
The model for the boiler unit of the plant was developed 
earlier I 1' i . In the present study an attempt is made 
to model the plant following the steam path do-vm-stream 
of superheater and back to drum through the regenerative 
feed heaters. The plant was also made available for open 
loop tests by the kind cooperation of the plant management. 
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1.1 LITERATURE REVIEW 

Investigations in the power station modelling 
started in late fifties when many investigators considered 
the problem of drum type boiler modelling. In 1958 
Chien et al,''_2] first presented a model and the dynamic 
analysis for a drum type boiler. The plant considered was 
an oil fired, single furnace natural circulation unit. The 
model consisting of a set of differential equations was 
solved using electronic differential analyser and a few 
open loop responses were obtained. But the model was too 
simple to represent the complicated dynamics of the boiler. 

Daniels ot al . '3 31 19^0 developed a model for 

the 200 MW unit of the Philadelphia Electric Co,, Crombay 
which was a piTLverised coal fired, twin furnace, controlled 
circulation, reheat boiler. Extensive open loop tests were 
carried out in the actual plant in order to validate the 
model. But the model did not produce sufficiently accurate 
responses to be considered suitable for control system 
design. 

In 1965, Thompson ' 43 derived a line '.rised mathe- 
matical model describing the dynamics of the drum, down- 
comer, waterwall, economiser, superi\(§ater and combustion sub- 
systems of a large utility power plant. Sufficiently 
accurate transient respons^-were obtained to permit the 
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perfomiance evaluation of proposed multivariable controllers. 
The model was verified by transient tests on the Crombaj?- 
unit of Philadelphia Electric Co, But this model was not 
developed for the full plant and was terminated at the 
throttle valve. Moreover, because cf large number of elemen- 
tal volumes, the model had a large state space dimension 
rendering it unsuitable for extensive control s3=^stem 
analysis, 

im report C 61 also presented a model for most 
general type of boiler and turbine using lumped parameter 
technique. Open and closed loop responses of the system 
were obtained and analog closed loop control operation 
was incorporated in the sim.ulntion model. Controller 
settings are determined using Ziegler - Nichols techniques 
based on single loop responses. 

In 1970, Mcdonald et al. [ 7 , 8 , 9 "j developed an 
extensive nonlinear mathematical model of a drum tjqpe, twin 
furnace, reheat, boiler- turbine-generator system. The 
results of an extensive programme of comparisions of closed 
loop steady state operation and open loop transient res- 
ponses of field data and the model are reported to 
substantiate the ■validity of the model. The relatively 
low order of the model, its high degree of accuracy in 
representing the relevant responses and its development 



7 


on the basis of a first principle analysis, make it ideal 
for control system analysis and design. The model para- 
meters are physical quantities v/iiich are readily obtainable 
from design data or unit acceptance test data. The model 
can be adapted to other applications also. 

In the present model the miiLtivalve approach 
followed by Mcdonald et al , i'9'! for the governing stage 
modelling is adopted. In the existing literature for the 
turbine modelling the medium pressure and low pressure 
turbines are lumped together and extractions are assumed 
to be taken at the end of the turbine. But, in the present 
model the medium pressure (MP) and low pressure (LP) tur- 
bines are modelled seperatcly. Moreover, the MP turbine 
is subdivided into four smaller turbines to accoimt for the 
reduced steam flow in the succeeding stages after each 
extraction. Furthermore, all the feedwater heaters are 
considered separately instead of Imping them into one 
feedwater heater as has been done upto now by other resear- 
chers. A brief outline of the content of each chapter is 
given below. 

In Chapter 2 plant description has been given. 
Flow path taken by main steam, condensate water and feed 
water has been presented and explained in detail. 



8 


In Chapter 3 the complete mathematical model of 
the plant has been developed. Further^ ass'umptions and 
simplifications made in the process of developing the model 
have also been listed. 

In Chapter 4 the results have been discussed. 

Plant responses have been obtained for a 10 M¥ load drop and 
responses to some variables have been compared with the 
results obtained from plant tests. Also steady state data 
from the plant has been compared with the steady state data 
of the model. F-iis chapter is concluded with a list of 
topics for research for further work. 



CHAPTER 2 


PLANT DESCRIPTION 

Uie Panki Thermal Power Station comprises of two 
old units of 32 MW each to which two nev/ units of 110 MI'J 
each have teen added recently. The second unit of 110 MW 
known as Unit No. 2 of the Panki Thermal Station Project 
is manufactured hy the Bharat Heavjr 3Illectricals Limited, 
India and was selected for modelling purposes. 

The unit consisting of a boiler and turbine is 
a single drum type, reheat, closed cycle system with rege- 
nerative feed heating. The boiler is a single drum, 
radiant water tube, dry bottom circulation, pulverised coal 
fired one with superheater and reheater. Since our study 
pertains to only the turbine part? the detailed specifica- 
tions are listed in Appendix 1 and tiie plant details are 
shown in Figure 1 . 

The flow of steam from the secondary supeih eater 
bifurcates into two equal halves and each half, after 
passing through a normally open throttle valve, enters 
a steam chest on either side of the turbine. Each steam 
chest contains four governing valves that are opened 
sequentially to control the flow of steam to the turbine 
as shown in Figure 2. Each valve discharges the steam 



10 


to its own section of the first stage nozzle ring. The 
valves are regifLated by secondary oil pressure which can 
be remotely controlled by a pump. This pump is driven by 
an electrical motor and can be controlled from the control 
room. Valves 1 and 2 open at the same secondary oil 
pressure. The different pressures ■ at which these 
valvesopen and are filLly open are given in Appendix 1. 

Two nozzle ring sections contain seven nozzles each and 
the other two eight nozzles each. The steam leaving the 
governing valve is expanded througli these convergent diver- 
gent nozzles and emerges at high velocity to impinge upon 
two rows of moving blades called the 'Curtis Stage' or the 
'governing stage'. This is followed by eight impulse stages 
of the high pressure turbine at the exit of which some 
steam is extracted for the high pressure feed water heater 
No. 2. The bulk of the exit steam from the H.P, turbine 
passes through a series of two reheaters called Trifflex 
and Exit reheater respectively. Exit reheater is in the 
first pass of the boiler and the T biflex in the second pass 
of the flue gases (not shown in figure). Steam from exit 
reheater with temperature controlled at 535 by attem- 
peration . is fed into the Mediutm Pressure Turbine which 
consists of twelve impulse stages. Steam is extracted at 
the end of 4th, 8th, 10th and I2th stages and is fed to 
High Pressure Feedwater Heater 1, Deaerator, Low Pressure 
Heater 5 and Low Pressure Heater No. 4 respectively in 



11 


that sequence. Steam from the exit of medium pressure 
turbine is fed directly into .a double-flow loxr pressure 
turbine which consists of four reaction stages on each side. 
Steam is extracted from both the sides at the end of 1st, 
2nd and 3rd stages and is fed into the low pressure feed- 
water heater 3} low pressure he. iter 2 and low pressure 
heater 1 respectively. The steam exiting from the lovr 
pressure turbine is condensed in two condensers each, and 
the condensate is passed through the series of low pressure 
feed heaters mentioned above. The steam condensed in the 
LP heaters 1 and 2 is fed back to main line before the feed 
water enters LP heater 1. Condensate from LP heaters 3? 
and 5 is fed into main floi\r lino at the entra^e of the 
LP heater 3* From the exit of the LP heater 5” the feed 
water is taken to the ’Feed Water Tanlc' through the deae- 
rator. From the feed water tank, water is taJ-on to the 
two high pressure heaters with the help of a food water 
pump. !jQie steam, condensed in the Wo heaters is fed back 
into the main line at the entrance of high pressure heater 
1. Exiting water from these heaters goes to the economiser 
where it is heated by the flue gases and is then passed on 
to the boiler drum. The details of flow path and other 
details from here onwards can be found in Reference [ 1 ’"' , 



aiAPTER 3 


FlATHMATIGAL MODELLING 

3,-] GEI^ERAL COMMENTS 

Panki Thermal Power Station details are given in 
three drawings [15 ^ • These being erection drawings they des- 
cribe the plant in great details wliich is not required for 
the present study. Hence, Figure 1 shows a simplified ver- 
sion of the turbine unit, the main steam path, feed water 
and condensate water lines wiiich x^as obtained from the above 
three referred drawings. Even "this diagram was too detailed 
because it included some features which will become opera- 
tional only in case of emergencies or turbine trips. Hence 
it was further simplified and Figure 3 shox/s the plant 
diagram as it was modelled. Steam admission to the high 
pressure turbine is controlled by four valves which admit 
steam to a nozzle ring followed by two rows of moving 
blades with an intervening row of fixed guide blades. These 
constitute the governing stage of tho turbine. The govern- 
ing stage is followed by eight impulse stages of the H.P. 
Turbine. After the H.P. Turbine some steam (approximately 
2 >%) is extracted for passing it through high pressure feed 
water heater No, 1 and the rest is passed through two 
reheaters called Triflex roheater and the Exit reheater 


as discussed in Chapter 2* 
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To obtain good representation for dynamic studies, the 
medium pressure turbine is subdivided into four smaller 
turbines with the l-umping made at the end of each extrac- 
tion. This results in first and second turbines having 
four stages each while third and fourth have two stages 
each. The extraction steam, from various extraction points 
in the M.P. Turbine, is taken to high pressure feed water 
heater No. 2, the feed water tanlc and low pressure feed water 
heaters No. 5 and 4 respectively in that order. Outlet 
steam from the medium pressure turbine? is fed into low 
pressure turbine. The arrangement of tliis turbine is sjmime- 
tric and as such, is replaced with a simple four stage 
turbine using the principle of symmetry. The steam that is 
being extracted at 3 points in the LP turbine is assumed to 
be extracted at the exit of each turbine section and is 
utilised for heating th.c feed water in heater Nos. 3} 2, and 
1 respectively. The two condensers, one each for the two 
sides of the low pressure turbine are also lumped into one. 
The condensate then passes through two low pressure feed 
water heaters which have been lumped into one as shown in 
Figure 3 , The bypass line to condensate storage tank 
is used only in emergency conditions and has not been 
modelled. The feed water then passes through feed water 
heater Nos. 2, 3 and k and then is fed into feed water 
tank. From this tank feed water is pumped into the high 
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pressure feed water heater No. 1 and then to the No. 2 
heater. From there the feed water goes to boiler drum 
through tile economiser. 

In developing the model, the system approadi has 
been followed and the components are considered to be 
blocks across which the physical changes in the condition 
of steam or water, as the case may be, are considered. 
These changes are represented in the form of mathematical 
equations by applying mass balance, momentum balance and 
energy balance across these blocks. The thermodynamic 
path of the process as it flows through the turbine is 
given in Figure 9. '^^le details of the model are presen- 

ted below. 

3.2 GOVEENING STAGE MODELLING 

There are four throttle valves to control the 
admission of steam to the high pressure turbine. The lift 
of these valves is controlled by the oil pressure. Out of 
these four valves, two throttle valves attain fully open 
position at the same oil pressure and hence these two have 
been lumped into one. As a result the steam admission has 
been modelled to take place through three valves which 
open sequentially as the load is increased. Following 
assumptions are also made: 



1 . steam condition at the inlet of the valve is same for 
all the three streams, 

2. steam flows in three streams caused by the three valves 
and undergoes different changes in steam condition, 
however pressure at the exit of governing stage is same 
for all the three streams, 

3. impulse blades ideally absorb all kinematic energy in 
the incoming steam and converts it into work. 


With the above assumptions, mass flow rate 
through the throttle va].ve is related to tlie pressure drop 
across it by the following equation [ 9 1 


2 





( 3 . 1 ) 


where 

(W^)y 




^SS 

^ss 


= flow through the throttle valve (Kg/ sec.), 

= flow area available for a partially open 
valve ( 

3 

= density of steam at the valve inlet (Kg/m ), 

2 

= steam pressure at the valve inlet (Kg/cm ), 

= nozzle ring pressure in three streams 

2 

respectively (Kg/ cm ), 

= a constants 
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Mass flow rate through the nozzle ring is given by 

= C p^j - Pj^p) (3.2) 


where, 

(Wp)]^ 

C 

^ THI 
%P 


flow through the nozzle ring (Kg/sec.), 
a constant, 

■3 

density of steam in the nozzle ring (Kg/m ), 
steam pressure at the exit of curtis stage 
(Kg/cm^) . 


Assuming no storage in the nozzle ring the f3.ows 
given by Eqns. 3»1 and 3*2 are equal, so we get the follow- 
ing equation for flow 


(¥^) 


2 


C. 




) 


(3.3) 


Here 


A^ C 




(Ay Pgg + Cy Pyjjy C) 


and Wn 


total steam flow at a particular 


For a fully open valve Ay = 1 and flow is 

similar equation 


load (Kg/sec.), 
given by a 


(Wy) 


*^2 ^SS ” %P^ 


where 


C 


THI 


( Pgg + Cy Pygy C) 


(3.^) 
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If the pressure ratio across the valve nozzle 
comhination is less than the critical <0‘So54) then 

sonic conditions exist and. flow is independent of the down 
stream conditions. In such a case the flow through a 
partially open valve is given [ 9 ] 

P 

= A (3.5) 

-L -b /rn 

/■^SS 

where, 

Tgg = temperature of steam at valve inlet 
= a constant. 


Flow through a filLly open valve is given (similar to Eq . 
3,5 with Ay = 1) by 


Wn 


^SS 


/'T 


( 3 . 6 ) 


ss 


Once the total flov7 through tho turbine, (W^,) is 
known, the preceeding equations are used to determine the 
pressure at the exit of the governing stage. But the 
situation becomes complicated if critical pressure ratio 
exists across the nozzle (i.e. ^°■5o51^for super- 

heated steam). In such cases any one combination of the 
following i.e. fully open subsonic, fully open sonic, 
partially open subsonic and partially open sonic conditions 
may occur. A generalised program has been developed to 
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take into account all such situations. The flo^; chart for 
this is given in Fig. ^ . First oil pressure is determined 
from the known load value (Fig, 5 ) and a va3-Ue of the 
governing stage exit pressure (Pgp) is assumed. The status 
of a particular valYe i.e. whether it is partially open or 
fully open is determined with the help of this oil pressure 
(Fig, 6 ). Now the partially open stream flow is deter- 

mined using Eqn. 3*3» For the fully open streams a check 
is made to determine whether critical pressure ratio 

exists across the valve nozzle combination. 
If critical pressure ratio does not exist, the flow is 
determined using Eqn. 3.4 while if critical pressure ratio 
exists Eqn, 3»8 is used. At this stage the sura of the 
three stream flows determined above is computed and is 
compared with the total turbine flow (VJp) , If the total 
of these individual flows is not equal to the total turbine 
flow (Wj), then the assumed P^p value is modified and the 
individual flows are determined again using this modified 
value. This procedure is repeated and the program is 
teiminated when we get the sum of three stream flows equal 
to the turbine flow within prescribed accuracy. 

At this stage a check is made to see whether 
all valves are fully open.. If the answer is yes, then 
the down stream pressure, temperature and enthalpy are 
determined using appropriate equations given below 
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(Eqns. 3*7 to 3*13). If all valves are not fully open, 
then a combination of partially open and fully open valves 
will exist. For the partially open valves the pressure 
^IHI nozzle ring is calculated using the subsonic 

flow equations 3.2. With this pressure known, the ratio 
is calcTilated. Hiis ratio then gives the number 
of sonic and subsonic partially open streams. Suppose 
subsonic flow exists in all the three streams then the 
above calculated sind P^p values are valid and the 

down stream pressure, temperature, enthalpy can be easily 
determined using Eqns. 3*7 fo 3.15* If a combination of 
sonic and subsonic flow streams exists then the subsonic 
stream flow is calculated with the values of Ppj.jj and Pp-p 
determined above while for the sonic streams is deter- 

mined by equating subsonic flow through the valve to sonic 
flow through its corresponding nozzle. Thus with this 


Pppi the flows through these sonic streams are also deter- 
mined using Eqn. 3.2. The sum of these sonic and subsonic 
flow is now evaluated and compared with the total turbine 


flow (Wp) . If these two flov/s are not equal, Ppp is again 
modified and the procedure is repeated until we get the 


sum of flows equal to the turbine flow within the desired 


accuracy. With the Ppp known the enthalpy and tempera- 
tures are determined using equations 3.7 to 3.1? descri- 


bed below. 
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Following are the state relations for determining nozzle 
ring steam temperatures obtained by using empirical fits 
from steam tables C161 


T 

where 

%1 

where 

%2 

^TH2 

where 

%H3 

^TH3 


0.5 + 1.33 - 620.0 (3.7) 

first stream nozzle ring temperature i°C) 
first stream nozzle ring enthalpy (Kcal/Kg). 

0.5 Pig2 ^'tFI2 " 

second stream nozzle ring temperature (°C), 
second stream nozzle ring enthalpy (.Kcal/Kg). 

0.5 Pq3J3 + 1.33 1%^ - 620.0 (3.9) 

third stream nozzle ring temperature (°C), 
third stream nozzle ring enthalpy (Kcal/Kg). 


Ideal exit temperatures of steam in the nozzles are deter- 
mined by using the isentropic laws given by the following 
expressions 


T 

■^HPI 



/ -"HP N n 

'' p ^ 

^THI 


( 3 . 10 ) 
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where 


first streaiH nozzle ring exit temperature ( C) 
index for super heated steam (n = 1,3). 


Prrp ^ 

p - T ( — ) 

■^HP2 “ I3i2 ^ 


(3.11) 


where 


Tj,jp 2 = second stream nozzle ring exit temperature (^C), 


T ( ^HP s 

"lT-13 P^i3 


(3.12) 


Nozzle ring exit enthalpy is determined using empirical 
fits obtained from steam tables for following state rela- 
tions . 


0.3 Pjjp + 0.617 + 527.6 


(3.13) 


where 


%P1 ~ enthalpy of steam for first stream at nozzle 

ring exit (Kcal/Kg), 


■0.3 PtTu + 0,617 T.fpo + 527.6 


(3.1^) 


where 


enthalpy of steam for second stream at 
nozzle exit, (Kcal/Kg). 


-0.3 Pjjp + 0.617 %p3 + 527.6 


(3.15) 
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where 

%P3 


enthalpy of steam for third stream at nozzle 
exit (Kcal/Kg). 


The average enthalpy of the three streams is determined by 
taking a weighted average of the enthalpies of the three 
streams using the energy balance equation. 


3.3 HIGH PKESSdKE TUBBIKE MODELLING 

Steam coming out of the governing stage enters 
the high pressure turbine. For modelling, it is assumed 
that high pressure turbine impulse blades ideally convert 
the kinetic energy into work. Furliier, it is assumed that 
the relationship between enthalpy, specific volume and 
pressure for the supeiheated steam, follows Callender's 
empirical relationship for all non~constant pressure 
processes [61. 

The flow through the high pressure turbine is 
related to the pressure drop across it by following equa- 
tion 110] 

(%) = K.J Pjjp (Pjjp - P(^j^) (3.16) 

where 

K.| = a constant to he detennined at initialisation 

time, 

2 

^CR ~ cold, reheat steam pressure (Kg/cm ), 
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The change in the outlet flow of turbine with a change in 
turbine inlet flow, as indicated by the following equation 
C 6 J , is determined by considering the turbine to have a 
first order transfer function. 


^d-XWcj,) 

di: 

where 


(1 - Cq) 


"0 



( 3 . 17 ) 


^CR ~ steam flow at the turbine outlet (cold reheat 
flow)(Kg/Sec.) , 

Cq = a constant, 

Tq = time constant for the high pressure turbine 
(Sec.) , 


Expansion of steam in the turbine is assumed to be 
isentropic (n =11,3) and the isentropic temperature of the 
cold reheat steam is given by the following equation Cl 01. 


'^CR' 

where 

^CR 

n 



p 

( CR 

' P 
^HP 


n--1 

n 


(3.18) 


temperature of the cold reheat steam (*^K), 

temperature of steam at governing stage 
exit (°K), 

index for the supeiheated steam. 


Value of 

noi. 


n 


is taken to be 1.3 for superheated steam 



Following state relation is used to determine the ideal 
cold reheat enthalpy in terms of cold reheat pressures 
and cold reheat temperature, 


= 0.593^ + 0.65 + 53^.069 

(3.19) 

where 

= cold reheat steam enthalpy (Kcal/Kg). 


The overall efficiency of the turbine is reduced 
by the effects of stage leakages, root and tip interference 
losses and rotational losses which increase witii the flow 
rate. These factors are grouped together in the model, 
in the form of a efficiency correction factor [ 6 ) 


n 


F 


(^) _ w 


R 


(^) (1 


W 


R 




C3.20) 


where 


lip = efficiency correction factor for the high 
pressure turbine, 

rated steam flow of the turbine (Kg/Sec.), 

= ratio of no load flow to the rated flow of 
the turbine. 

Variation of efficiency correction factor, with 


the ratio of actual turbine flow to the rated flow, is 
shown in Fig, 8, 
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The overall efficiency is given by 


^HPA “ ’^F 


( 3 . 21 ) 


where 


Tlg.p = thermodynamic efficiency. 
The actual reheat enthalpy 


(Referring to Fig. 9 ) 


is given by 


^CEA. = %P ■ "'hPA ^ %P " 


( 3 . 22 ) 


where 


= actual cold rehea 


it stoam enthalpy (Kcal/Kg) 


The power produced by the HP turbine in terms 
of the flow rate and enthalpy drop in the HP turbine is 


given by 


POVJ^ = 0.0042 ^MPA ^%P ” (3*23) 


where 


POW^ = power produced by the high pressure turbine 
(MW). 


3 .If REHEATER MODELLING 

The following assumptions are made: 

1. The average values of the parameters like enthalpy, 
temperature and density are assumed to be the exit 

values. 

2. The steam leaving the Irlflex reheater enters the 
Etlt reheater and has the instantaneous properties 
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of the steam leaving the Triflex reheater. 

3. Steam is heated mainly by convection. 

4-, Reheater tubes are replaced by eqirivalent control 
volume across which mass balance energy balance and 
momentimi balance are considered, 

5. Heat transfer rate from the flue gases to the tubes is 
considered constant. This is actually not so, however 
if the boiler dynamics is also included then this will 
be determined by that dynamics. 


3.4-.1 Triflex Reheater 


The following equations arc used. 

Mass Balance 

Mass balance across the Triflex rcheatcr is given 
by the following equation [791 


^^CR ” ^^R01 


■ d ( P 


dt 


^Rl) 


( 3 .^) 


where 

Wroi = Triflex rdaeater exit steam flow rate (Kg/Sec.) 
p _ . = density of the steam at the Triflex reheater 

xtT 

o 

exit (Kg/m'^), 

"^Rl " steam volume of the Triflex reheator (m ) . 
Momentum Balance 

The monentum equation describing the steam 


pressure drop is 
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^CR “ ^E01 


2 

R1 Pr^ g 


d_^) 

dt 


(3.25) 


where 



= friction factor for the Triflex reheater, 
= effective length of the pipe (m), 

= area of the pipe (m ) . 


Substituting the differential tem in equation 3.25 from 
equation 3.17 we get 


V/ 


^CR " ^ROI 


R1 


CR 

) 

R1 


Xj . 


(1 - Cq) 


A. 


T, 


0 


( 3 . 26 ) 


Energy Balance (Gas Side) 

The energy storage in the metal tiibo walls is 
given by the following equation [1 3, 


‘^RGI “ *^RS1 ~ ^ 


IM1 


M. 


dt 


(3.27) 


where 


^RGI 

‘^RSI 

'^IMI 

^RM1 


= heat flux on the tube walls in Triflex 
reheater (Kcal/Soc.), 

= heat added to the steam in the Triflex 
reheater (Kcal/Sec.), 

= specific heat of the Triflex reheater 
tube metal (Kcal/Kg °C), 

= metal mass of the Triflex reheater (Kg), 
= tube metal temperature in the Triflex 
reheater (°C). 
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Energy Balance (Steam Side) 


The energy balance for steam in the Triflex 
relieater is given by the following equation 


1 

Y 


dt 


^E01^ 


where 


“^RSI " ^^B01 ^E01 ^■^CR ’^CR 

( 3 . 28 ) 


M 


R1 


effective st -eam mass in the Triflex rehcater 


(Kg), 

Y = ratio of specific hearts (Cp/Cy for steam). 
Following state relation is determined in terms of Triflox 
reheater exit pressure and enthalpy using emplri' 

cal fits from the steam tables [16] . 


Tpgi = 0.78? (PpQ^ - 2h,h7) + 1.81 (Hp^^ - 772.0) + 397 


(3.29) 


where 


^RSI ~ temperature of steam in the Triflex reheater ( C; 


The heat added to the steam from the metal walls is given by 
the following equation C 1 ] 


^RS1 “ ^EMS1 ^^RM1 ” ^RS1^ 


(3.30) 


where 


^MS1 ~ ^ constant to be detcimined at the initia- 


lisation time 



29 


Because of the large volume involved in the 
reheater there is a storage efiect in it. As a result we 
can consider it as a tank across which changes in flow are 
approximated hy a first order transfer function. The 
following equation is therefore used 


d 

dt 




VJ - ■'■•j 

'^512 

k' 


( 3 . 31 ) 


where 

Ty = time constant associated vrith the exit 
reheater (Sec.), 


3.^.2 Exit Rehoater 


The assumptions and arguments made in the case of 
Triflex reheater. Art. are valid in this case also 

and hence will not he restated. 


Mass Balance 

Mass balance across the exit reheater is given by 


^EI2 ^ ^ED2 


dt 




( 3 . 32 ) 


where 


^EI2 

^ED2 

^E2 


= steam flow entering the Exit reheater (Kg/Sec.), 

= steam flow at the exit of Exit reheater 
(Kg/Sec.), 

= density of steam at Exit reheater exit 
(Kg/m^) , 
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3 

^R2 ~ steam volume In the Exit roheater ( m ) . 

Momentum Balance 

Momentum balance describing the steam pressure 
drop in the Exit rehoater is given by the following 
equation 

p p _ p , ^2 d_ ^ 

^EI2 " BD2 “ R2 Pj^2 g dt '^‘■^RI2^ 


where 


■EI2 


ED2 


F 


W. 


R2 

KC2 


pressure of inlet steam to the Exit 
reheator (Kg/cm^), 

pressure of steam at the exit of Exit 
rehoater (Kg/cm~), 

friction factor for the Exit rehoater, 
irfLet flow rate oi' steam to the Exit 
reheatcr (Kg/sec.), 

L 2 = effective length of pipe (m), 

2 

A 2 = area of the pipe (m ). 

Energy Balance (Gas Side) 

The energy storage in the metal tube walls is 
given by the following equation 


*^^12 “ ^RS2 ~ ^TM2 ^2 d?f ^^RM2^ 


(3.34) 


where 


^EG2 


heat flux on the tube walls in Exit 
reheater (Kcal/sec.), 



^RS2 

^m2 

M2 

^iM2 
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heat added to the steam In the Exit reheater 
(Kcal/ sec, ) , 

specific heat of the Exit reheater tube 
metal (Kcal/Kg °C), 

metal mass of the Exit reheater (Kg), 
tube metal temperature in the Exit reheater 

(°c). 


Energy Balance (Steam Sid e) 

Energy balance for steam is given by the follow- 
ing equation 

^ ^ ^R2 ^%02^ “ %S2 ^'EI2 ^RI2 ^^E02 ^BD2 

(3.35) 

v/here 

^R2 ~ effective steam mass in the Exj.t reheator (Kg), 

^HE2 ~ enthalpy of steam entering the Ejcit reheater 
(Kcal/Kg), 

^BD2 ~ enthalpy of steam coming out of the Exit 
reheater (Kcal/Kg). 

Following state relation is used in determining Exit 
reheater temperature ^sing empirical fits from 

steam tables fl6] . 

^HS2 = ^^B02 “ 23.12) + 2.05 - 8^+5. 0) 

+ 535.0 ^ (3.35A) 
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Ihe heat added to steam in the Exit reheater is given by 
the following expression fl ] 


Q 


ES2 


= C 


IMS2 




0.8 


(1 


IM2 


^ES2^ 


(3.37) 


where 

*^11102 ~ ^ constant to be determined at the time of 

initial! zation . 


Because of the big reheator volume there is sto- 
rage effect in the reheater. So we can consider the 
reheater as a tank across which changes in f3.ow are approxi 
mated by a first order transfer function. The following 
equation is therefore used 


dt 




^RI2 " ^^B32 
■^8 


(3.38) 


where 

Tg = time constant associated with the Exit reheater 
(Sec.) . 


3.5 MEDIUM PRESSURE TURBINE MODELLING 

The medium pressure turbine is considered to 
consist of four smaller turbines to account for the reduc- 
tion in flow rate to the succeeding Impulse stages on 
account of extractions at the end of these stages. 
Extractions are now assumed to take place at the end of 
each turbine. In addition, inlet flow to the M.P. turbine 
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is supplemented by attemperation flow to control the 
temperature* Ihe condition of steam to the HP turbine is 
therefore modified by the attemperation flow. This is 
taken care of in the model by using mass balance and energy 
balance to the attemperator to give the flow conditions at 

MP turbine inlet. 

Inlet flow to the Medium Pressure (M.P.) turbine 
is given by the fallowing equations 


= ^m2 %iMP (3.39) 

where 

Wj^pi = inlet flow to the M.P, tiirbine (Kg/sec,), 

^^ATMP ~ attemperation flow required to bring the 

steam to a specified temperature (Kg/sec), 


Energy balance for attemperator is given by the expression 


^ED2 ^^B02 “ %PI ^MPI \tMP %TMP 


where 


%PI - 
%TMP = 


enthalpy of steam entering the M.P. turbine 
(Kcal/Kg), 

enthalpy of the attemperation water 
(Kcal/Kg). 


Once the flow rate to th.e M.P, turbine is known 
the down stream conditions can be determined using the 
same equations as were used in H.P. turbine (Art, 3*3), 



thus the flow through each turbine is related to their 
respective inlet and outlet pressures by following equa- 
tions [ 9 ]• 


3 ^ 


2 

%PI = ^2 ^ ^ %P1^ 


(3.^1) 


where 


%PI 

%P1 


= ' a constant to be determined at the time of 
initialisation, 

= pressure of steam at the inlet of first M,P, 

p 

turbine (Kg/cm ), 

= pressure of steam at the exit of first M.P, 

2 

turbine (Kg/ cm ) . 



^3 ^MP1 ^%P1 


- 


(3/^2) 


where 

^3 

^MP1 

%P2 

PMP.1 


a constant to be determined at the time of ini- 
tialisation, 

stean flow rate entering the second M.P, 
turbine (Kg/sec,), 

pressure of steam at the exit of second M.P, 

2 

turbine (Kg/cm ,), 

density of steam at the exit of first M.P, 
turbine (Kg/m"^), 
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V/, 


h^P2 ~ % PmP2 


%P3^ 


(3.^3) 


■where 


= 

%P2 ~ 

%P3 = 

P 

MP2 “ 


a constant to be determined at the time of 
initialisationj 

steam flovz rate entering the third M.P, 
turbine (Kg/sec . ) , 

pressure of steam at the exit of third M.P. 
turbine (Kg/cm ) , 

density of s'team a'b the exit of second M.P, 

O 

turbine (Kg/m"^). 



^5 ‘^MP3 ^%P3 


- %po^ 


(3M) 


where 


% 

%P3 

%po 

^MP3 


= a constant to be determined at the time of ini- 
tialisation, 

= steam flow rate entering the fourth MP turbine 
(Kg/sec.), 

= pressure of steam at the exit of four-thM.P. 
turbine (Kg/cm ), 

= density of steam at the exit of ■third M.P. 

O 

turbine (Kg/m"^). 

The change in the ex:it flows as a first order 


approximation of a change in each turbine inlet flow is 
described by the foliowing equations i 
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dt 




where 


(1 - C^) - w, 


MP1 


(3 > 5 ) 


= a constant for deteimining extraction flow, 

= time constant of the first M.P, turbine (sec.), 


^ (IT _ (1 C^) 

dt '^"mP2^ ~ T3 

where 


( 3 .^ 6 ) 


= a constant for determining extraction flow, 

T3 = time constant of the second M.P, turbine (sec.). 


dt 


(W^P3) 



" *^6^ ''SlPP “ 
\ 


( 3 .^ 7 ) 


^diere 

= a constant for determining extraction flow, 

Ti_|^ = time constant of the third M.P. turbine (Sec,), 


^ (u ) - (1 - Gy) W^p3 -> W^pQ 

dt ^^MPO^ - 

where 


( 3 . 4 - 8 ) 


Cy = a constant for determining extraction flow, 

= time constant of the fourth M.P. turbine (sec.). 

Expansion of steam in the four turbines is 
assumed to be isentropic and the following equations are 
used to determine the isentropic temperatures. 
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T _ rp / MP1 S 
-^MPI “ -^MPI Pj^pj ^ 


n-1 

n 


( 3 .^ 9 ) 


where 


%PI ~ temperature of inlet steam to first M.P. 
turbine (°K), 

%P1 ~ temperature of steam at the exit of first 

M.P. turbine (°K) , 


^>IP2 ~ ^MP1- ' P,,, 


n-j 

(^7 


MP1 


( 3 . 50 ) 


where 


%P2 ” temperature of stecun at the exit of second 


M.P. turbine (*^K) . 
^MP3 - ^MP2 ^ Pj^pp ^ 


n-1 

n 


( 3 . 51 ) 


where 


Tj^P 2 = temperature of steam at the exit of third 
M.P. turbine (‘^K). 


%P0 


^MP3 


n-1 

^MPO 

%P3 


( 3 . 52 ) 


where 


m 

■^MPO 


temperature of steam at the exit of fourth 
M.P. turbine (°lO. 


The following empirical fits are used to get 
the enthalpy of steam at the exit of each M.P. turbine 
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%4P1 “ ~o»35 0*52 + 571 tS (3*53) 

■where 

%IP1 ~ enthalpy of steam at the exit of first M,P, 
t nrbi ne ( Kcal/Kg ) . 

^MP2 ~ ••0,5 0»51 -^p[p2 58o ,4 (3»54) 

"where 

%P2 ~ enthalpy of steam at the exit of second M.P. 

turbine (Kcal/Kg), 


H 


'MP3 ■” -0.636 2j,^p 2 0.5 T'^jp^ 584,5 (3*55) 


■where 


%1P3 ~ enthalpy of steam at the exit of third M.P. 

turbine ( Kcal/Kg ) . 


H, 


where 


B 


■MPO 


■MPO = ’•'•*36 P^po + 0,48 + 592.18 (3.56) 


enthalpy of steam at the exit of fourth M.P, 


turbine (Kcal/Kg). 

Efficiency correction factors for the four turbines are 
deteimined as follows i 


El W, 


/ %P1n „ 

(-^) . Wp 


where 


( MP1 w. 

( ^ )(1 - Wp) 


(3.57) 


:= efficiency correction factor for first M.P, 
turbine, 
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ts = 


- w, 

- »L> 


(3.58) 


where 


"i2 = 


efficiency correction factor for second M.P, 
turbine . 


"§3 = 


•H 




(3.59) 


where 


F3 


efficiency correction factor for third M.P. 
turbine . 


/^MPO^ „ 


'F4 - H 


/ MPO \ / ^ T T \ 

- ''l> 


( 3 . 60 ) 


where 


efficiency correction factor for fourth M.P, 
turbine . 


The actual efficiencies are determined by the following 
expressions 

%P1 == ’^MP ^F1 (3.61) 

\fdiere 

= actual efficiency of the first M.P, turbine, 

=: theimodynamic efficiency of each M.P, turbine. 



4-0 


’^MP2 ~ 


( 3 . 62 ) 


where 


^MP2 ~ actual efficiency of the second MP turbine, 

%P3 = ’^MP ’^F3 (3.63) 

where 

= actual efficiency of the third M.P. turbine, 
^MP4- ~ ’^MP (3*6^) 


where 


^MPO ~ actual efficioncy of the fourth M.P, turbine. 

Actual enthalpy of steam is determined using the following 
equations . 

%P1A = %PI *■ ’^MPI ^^MPI ” ^P1^ (3.65) 


where 


%P1A ~ actual enthalpy of steam at the first M.P, 
turbine exit (Kcal/Kg), 

%P2A = %P1A ■“ ^MP2 ^^SdPlA ” %P2^ (3.66) 


where 


®MP2A ~ actual enthalpy of steam at the exit of 
second M.P, turbine (Kcal/Kg), 

^MP3A %P2A ” ^MP3 ^%P2A " %P3^ (3.6?) 


where 


= actual enthalpy of steam at the exit of 
third M.P. turbine (Kcal/Kg). 



where 


%POA 


%P3A “ ^MPO ^^MP3A “ %P0^ 


( 3 - 68 ) 


%POA ~ actual enthalpy of steam at the exit of 
fourth M.P. turbine (Rcal/Kg). 

Power generated by each turbine is given by the following 
expressions s 

POM 2 = 0.0042 (3.69) 

where 

P0¥2 = power produced by the first M.P. turbine (MW), 

POW^ = 0.00^2 n^jp^ ‘‘^MPIA ” %P2^ (3-70) 

where 

POW^ = power produced by the second M.P. turbine (MV7) 

P 0 W]_^ = 0.004-2 Wj^pp n^^p^ ^^P 2 A '■ %P 3 ^ ( 3 . 71 ) 

where 

POW)^ = power produced by the *tbird M.P, turbine (MW), 
POWj = 0.0042 w^3 "MP4 %P3A - (3.72) 

where 

POW^ = power produced by the fourth M.P, turbine (MW) 

3.6 LOW PRESSURE TURBINE MODELLING 

Exit steam from this turbine goes to the conden- 
ser. Extractions are lumped into one extraction and this 
extraction is assumed to take place at the end of this 
turbine. Further condenser pressure is ass-umed to be 
constant. 



The change in the Low Pressure (L.P.) turbine 
exit flow is assumed as a first order approximation of the 
inlet flow and is given by the follotring equation 

d 


dt ^^LPO^ ~ 


(1 - Cg) - WppQ 


"6 


( 3 . 73 ^ 


where 


Uppo = outlet floij from the L.P. turbine (Kg/sec.), 
Tg = time constant of the L.P. turbine (sec.). 


In addition efficiency correction factor, enthalpy of steam 
at the exit of L.P. turbine and power produced by the L«P. 
turbine are given by tlie following expressions. 

) - Wp 

(3.7^) 

) (1 “ VJj^) 

where 

= efficiency correction factor for the L.P, 
turbine . 


w. 


LPO 




¥ 


R 


W. 


(H 


LPO 


¥■ 


R 


’^LPA = ^LP 
where 


(3.75) 


\PA 


T1 


LP 


actual efficiency of the L.P. turbine^ 
thermodynamic efficiency of the L.P, turbine. 


%POA ~ ^MPOA ~ ^LPA ^%POA "" ^LPO^ 


(3.76) 



^3 


where 


^LPOA “ 
®LPO “ 


actual enthalpy of steam at the exit of 
L.P. turbine (Kcal/Kg), 

Ideal enthalpy of steam at the exit of L.P. 
tu rbine ( Kcal/ Kg ) , 


Power produced is given by the following expression, 


POWg _ 0.0042 npp_^ *^®MP0A “ ^LPO^ ( 3 . 77 ) 

where 

POW^ = power produced 'bj the L.P. turbine (MW). 

3.7 CONDENSER MODELLING 

It is assumed Wiat condenser pressure remains 
constant. Heat rejected by the condensing steam is there- 
fore determined by the quantity of condensing steam. 

Heat flow rate from metal to cooling water is 
given by the following expression [ 1 ] . 

%W = ^CMW “ ^CW^ (3.78) 

= heat flow rate from metal to cooling water 
in condenser (Kcal/sec,), 

= a constant determined at initialisation time, 

= quantity of cooling water (Kg/sec.), 

= ten^erature of condenser metal (°C), 

= tanperature of the cooling water (^^C), 


^ere 


%W 


'CMW 


W 


cw 


T 


'CM 

CW 



Enthalpy at the condenser exit Is determined hy an energy 
balance equation and is given by the following expression 


^co “ ^Lpo / '■'^LPO 


(3.79) 


where 


I'co 

Qsm 


= enthalpy of water at tlie condenser exit 
(Kcal/Ivg), 

= heat rejected by steam to the cooling water. 


Energy storage in the metal walls is given by the follomng 
equation 

^GM ^3 ^ (3.80) 

where 

“ specific heat of the condensor metal (Kcal/Kg *^0) 
= condenser metal mass (Kg). 


3.8 LOW PEESSUEE (L.P.) AKD HIGH PEESSUEE (H.P.) 

HEA-TEE MODELLING 

The feed water heaters are considered seperately 
to take into account the fact that extraction rates are 
different and they have different heat contents. For each 
feed water heater, energy storage in the metal tube walls, 
heat transfer rate from metal tube bo feed water are con- 
sidered and thus the enthalpy of feed water coming out of 
the heater is determined. 



3.8,1 Low Pressure Feed Water Heater No, 1 

Energy storage In the metal walls of the feed 
water heater No. 1 is given by the following equation 


%I1SM ” \hiMW 


^HMI \ dt 


(3.81) 


where 


'^LHiai 


‘^LHIMW 


%M1 

Mi^ 

^HM1 


= heat transfer rate from the steam to 

metal in the first L.P, feed water heater 
(Kcal/sec . ) , 

= heat transfer rate from metal to feed water 
in the first L.P, feed water heater 
(Kcal/sec .) 5 

=r specific heat of the tube metal in the 
first L.P. feed water heater (Kcal/Kg °C), 

= first L.P, feed water heater metel mass (Kg), 
= first L.P, feed water heater metal tempe- 
rature (°C), 


Hetal to feed water heat transfer rate in the 

heater is given by the following equation 1 11. 

0,6 


^LHIMW 


1MW1 




(T. 


mil 


T ) 

mwr 


(3.82) 


where 

%MW1 

^LHI 


a constant to be determined at initialisa- 
tion time, 

feed water flow rate through 'bhe first 
L.P, feed water hoator (Kg/sec, ), 
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^HWI = feed water temperature in the first L.P* 
feed water heater (°C) 

Enthalpy of feed water coming out of the first feed water 
heater is deteraiined by 

^LHIO ~ ®C0 ^ (3.83) 

where 

^LH10 “ enthalpy of feed water coming, out of the 

first feed water heater (Kcal/Kg) , 

3.8,2 Low Pressure Feed Water Heater No, 2 Modelling 

Energy storage in the tube metal walls is given 
by following equation 

^LH2M¥ “ *^2 * (3.84) 

= heat transfer rate from steam to metal in 
second L.P, feed water heater (Kcal/Sec.), 

= heat transfer rate from metal to water in 
the second L.P. feed water heater (Kcal/Sec.) 
= specific heat of the tube metal in second 
L.P. feed water heater (Kcal/Kg °C), 

= second L.P, feed water heater metal mass 
(Kg.) 

= second L.P. feed water heater metal tempe- 
rature (°C). 


^LH2SM 

where 

*^LH2SM 

%H2MW 

%M2 
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Metal to feed water heat transfer rate in the heater is 
given by following expression. Ll_ 


‘^LH2M¥ 


^HMW2 * • ^^HM2 ~ %W2^ (3.85') 


where 

= S' constant to be determined at initialisa- 
tion time, 

^LH2 ~ feed X'/ater flow rate through the second L.P, 
feed water heater (Kg/Sec.) 

%W2 “ feed water temperature in the second L.P, 

feed water heater (*^0) 


Enthalpy of feed water coming out of the second feed water 
heater is determined by 


%H20 “ ^LHIO ^ '^LH2MW ^ ’'^LH2 (3.06) 

where 

Hlh 2 o “ enthalpy of feed water coming out of the 
second feed water heater (Kcal/Kg) 

3.8.3 Low Pressure Feed Water Heater No. 3 Modelling 

Energy storage in the tube metal walls is given 
by following expression 

*^LH3aM ” ‘^ni3MW ~ ^3 • ^6 • ft (3.87) 

where 

Qlh 3 sm = heat transfer rate from steam to metal in 

the third'L.P. feed water heater (Kcal/Sec.), 



^LH3MW 

^6 
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= heat transfer rate from metal to water in 
the third L.P. feed mter heater (Kcal/Sec.), 
= specific heat of the tube metal in the third 
L.P. feed water heater (Kcal/Kg °C) 

= third L.P, feed water heater metal mass (Kg.) 
= metal temperature in the third L.P. 

feed water heater (°C) 


Metal to feed water heat transfer rate in the third heater 
is given by following expression 

0.6 

^LH3MW ~ ^HMVG '* ^^LH3^ • ^%13 " (3.88) 


where 

*^HMW3 

^LH3 

T 

HW3 


a constant to be determined at initialisation 
time, 

feed water flow rate through the third L.P. 
feed water heater (Kg/Sec.) 

feed water temperature in the third L.P, feed 
water heater (°C) 


Ehthalpy of feed water coming out of third feed water 
heater is deteimined by the following expression 


^LH30 “ ^LH20 *^LH3MW ^ %I3 

where 

^LH30 ^ 


(3.89) 


enthalpy of feed water coming out of the 
third L.P, feed water heater (Kcal/Kg). 



3.8.*+ Low Pressure Feed Water Heater No. ^ Modelling 


Energy storage in the tuhe metal walls is given 
hy the following expression 


■■ *^LH4MW 


^mh * ^^7 


dt 



) 


(3.90) 


where 


%Hi+SM 


= heat transfer rate from steam to metal in 

the fourth L.P, feed water heater (Kcal/Sec.), 
= heat transfer rate from metal to water in 

the fourth L.P, feed water heater (Kcal/Sec.), 
= specific heat of the 'tube metal in the 

fourth L.P, feed water heater (Kcal/Kg °C), 

= fourth L.P, feed water heater metal mass (Kg), 
■= metal temperature in the .fourth L.P, feed 
water heater (*^0) 


Metal to feed water, heat transfer rate in the fourth heater 
is given by .the following expression ClU 






0,6 


(T, 


■HM4 “ hiW4 


) (3.91) 


where 


%M¥4 

Wlh4 

%W4 


a constant to be determined at initialisa' 
tion time, 

= feed water flow rate through the fourth 
L.P. feed water heater (Kg/Sec.), 

= feed water temperature in the fourth L.P, 
feed water heater ('^C). 
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Enthalpy of feed water coming out of the fourth heater is 
determined by the following expression 


^LH^O 


where 


^LH40 


^LH30 ^ (3.92) 

enthalpy of feed water coming out of the 
fourth L.P. feed water heater (Kcal/Kg). 


3 . 8.5 High Pressure Feed Water Heater No. 1 Modelling 


Heat energy storage in the metal tube wall is 
given by the following expression 

d 


Q 


‘HH1SM 


%IH1M¥ “ %15 * % * dt ( 3 . 93 ) 


where 


Q- 


'HH1SM 


%H2MW “ 


%M5 


heat transfer rate from steam to metal in the 
first H.P, feed water heater (Kcal/Sec, ), 
heat transfer rate from metal to water in the 
first H.P. feed water heater (Kcal/Sec.), 
specific heat of the tube metal in the 
first H.P. feed water heater (Kcal/Kg ^C) , 
first H.P, feed water heater metal mass 
(Kg.), 

first H.P, feed water heater metal tempe- 
rature (°G). 

Metal to feed water, heat transfer rate in the first H.P, 
heater is given by the following equation Til. 


M 


8 


^HM5 
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%H1MW 


%1W5 * 


0,6 


• " ^HW5^ 


where 


"HMW5 




HH1 


T- 


HW5 


a constant to be determined at initialisa- 
tion time, 

feed water flow rate through the first H.P. 
feed water heater (Kg/Sec.), 

feed water temperature in the first H.P. feed 
water heater (°C), 


Enthalpy of feed water coming out of the first H.P. heater 
is determined by the following expression 


%I10 “ \eU-0 %I1M¥ ^ %I1 


(3.95) 


where 


^■HHIO “ enthalpy of feed water coming out of the 
first H.P, feed water heater (Kcal/Kg), 

3*8,6 High Pressure Feed Water Heater No. 2 Modelling 

Heat energy storage in the metal walls is given 
by the following expression 


Q. 




^HH2MW “ ^6 * ^9 • dt ^^HM6^ 


(3.96) 


where 


%iH2SM “ heat transfer rate from steam to metal 

in the second H.P, feed water heater (Kcal/ 

Sec.), I.I.f. KANPUR 

Hu. ■. rL 
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%H2MW 

^6 

M 

g 

%M6 


= heat transfer rate from metal to water in 

the second H.P, feed water heater (Kcal/Sec.), 
= specific heat of the tube metal in the second 
H.P. feed water heater (Kcal/Kg °C), 

= second H.P, feed water heater metal mass (Kg), 
= second H.P. feed water heater metal tempe- 
rature (°C). 


Metal to feed water heat transfer rate in the second H.P, 
heater is given by the following equation ■ C1J- 

^HH2MW ^HM¥6 * * ^%M6 " ^HW6^ (3.97) 

where 

= a constant to be determined at initialisation 
time, 

%H2 ~ feed water flow rate through the second H.P. 

feed water heater, 

%W6 = feed water temperature in the second H.P, 

feed water heater (^C)^ 

Enthalpy of feed water coming out of the second H.P feed 
water heater is given by the following equation 


®'HH20 

where 

%ffl20 


^HH10 ^HH2MW ^ %I2 


(3.98) 


enthalpy of feed water coming out of the 
second feed water heater (Kcal/Kg). 



3.8.7 Model Structure 
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The structure of non-linear ordinary differential 
equations is of the type 

X = f ( X, y, u ) (3.39 ) 

where x is the state vector given hy 
T 



^R1’ 

®R01’ 


^EI2’ ^'R2’ ^^E02» 

*'^E02’ 

^IM2’ 

%P1 


^''MP3» 


^HM2’ 

^HM3» 

^HM4- 

j 





y the auxilliary variables are given by 
y = G ( X, u ) (3.100) 

U is the control vector defined by 

U = (W^, Qp^Q2 ^ 

Substituting for T from eqixation ( 3 . 100 ) in 
equation (3«99) "we get the equation (3*99) as 

X = F < X, U ) ' (3.101) 

where F denotes the new functional dependence. 

For the above set of equations a computer program 
was developed and these equations were solved, on a digital 
computer (IBM 704-4) using the Piinga-Kutta integration < 
technique. The transient re.sponses of the turbine sys- 
tem for step changes in few control parameter were obtained 
and have been discussed in the next chapter. 



CONCLUSIONS 


In this chapter the dyinamic model of the turbine 
unit using physical laws and empirical relations was deve- 
loped. It yielded 22 first order differential equations 
alongwith 77 algebraic equations to represent the turbine 
unit. The above system of equations was solved in a digital 
computer and responses obtained will be presented in 
Chapter 4 -. 



CHAPTER if 


RESULTS 


lf.1 INITIALISATION 


Mathematical model for the Pajiki turbine was 
presented in Chapter 3. Before any transient study is 
conducted on the model the steady state of the model has 

to be obtained for given input conditions. This is called 
as initialisation. 


The constants used in the equations are ""determined 
using the steady state values of parameters at 90 MW and 
are given in Appendix 1 . For any given input condition 
there are two methods for determining the steady state of 
the model. In the first method the X • teim of equation 
3.99 is equated to zero and the resulting nonlinear alge- 
braic equations are solved by any suitable numerical 
elimination technique. In the second method, the steady 
state data obtained from the plant is used as a starting 
point in the integration scheme. The second method is 
used in the present study. As mentioned earlierj steady 
state data at 90 MW was available from the plant and this 
was used to determine the model constants. It can be 
pointed out here that generally the actual steady state 
data obtained from the plant may not be the same as that 
given by the model. However, in the present case, as 
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the steady state constants of the model were detemined 
from the plant data, the model steady state is same as the 
plant steady state. Fourth order Runga-Kutta method is 
used for the integration and the integration scheme 
is tenniiiated \j-hen the difference in the successive values 
of the state variables is reduced to less ttian of the 
previous values. Steady state operating values of the 
model are determined for many loa-ds and are presented for 
the case of the governing stage. In one case, for v/hich 
the plant data was available, comparison has been made with 
the plant data. Transient tests are carried out for two 
cases. First, for the 10 MW load drop from ^0 MW level 
and second for a step change in gas to metal heat transfer 
rates in the Triflex and Exit reheater. 

4-. 2 STEADY STATE RESULTS 

Table 4.1 shows the model steady state values 
obtained for 70, 80, 90, 100 and 110 !1W loads. The para- 
meters shown are pressures and temperatures after the three 

throttle valves P^p* ^IH3» %I11’ ^TR2’ 

obtained for the governing stage after implementing the 
flow chart shown in Figure 4, As tine load increases tine 
pressure after the three throttle valvesalso increases. 

This is explained by the fact that with an increase in 
load the throttle valve opening also starts becoming 
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larger and larger and as a result the pressure drop, across 
the throttle valve decreases. This has been shown in 
Table h,2. Further, for any 'one particular load, the pre- 
ssure after the first throttle valve will have a higher 
value than the second and the second will have a higher 
value than for the third. This is so because the first 
valve starts opening the earliest and as such at any par- 
ticular load it will have a larger opening than the second 
and similarly the second valve will have a larger opening 
than the third. The temperatures after the tiirottle valves 
also show the same trend and are presented in Table 4,1. 

TABLE 4.1 

PREdSUEES AITO TEMPERATURES AFTER THROTTLE VALVE 
FOR DIFFERENT LOADS 


LOAD 

PARAMETER 

! 1 

1 70 ! 

f r 

t 

80 i 

f 

90 I 

( 



p 

TH1. 

78.27 

92.66 

102.05 

118.79 

130.0 

^TH2 

68 .78 

8 1 #4*0 

88.55 

109,99 

118.66 

^TH3 

NOT OPEN 

NOT OPEN 

NOT OPEN 

107.19 

115.49 

T 

•^THI 

51^.00 

520.70 

5^.50 

531.20 

535.00 

^TH2 

511.20 ' 

516.10 

519*10 

527.60 

531.00 

^TH3 

NOT OPEN 

NOT OPEN 

NOT OPEN 526.50 

529.80 
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TABLE 4-. 2 

DBOP ACROSS THE "^AL^'ES from CONSTANT INLET 
PEESSUEE OF I 30 KG/CM^ 


LOAD 

PARAMETER 

(MW) j 

I 

! f 

i 80 i 

t ] 

! 

90 i 

t 

0 

0 

110 

^TS1 

51.73 

37.34- 

22.0 

11.21 

0 

^TH2 

61.20 

!+9,00 

4-1 .5 

20,00 

1 1 .4- 

^TH3 

NOT OPEN 

NOT OPEN 

NOT OPEtT 

23.00 

15.0 


The governing stage pressure %P and temperature 
Tjjp have heen deteimined for different loads and are pre- 
sented in Table 4-,3» The steady state value of P^p at 
various loads has been compared with plant steady state 
value at these loads in the same table. 


TABLE 4-. 3 

GOVERNING STAGE PRESSURE AT VARIOUS LOADS 


PARAMETER 

LOAD (MW) 

! 

f 

70 i 

1 

f 

80 { 

r 

r 

90 i 

1 

t 

100 i 

f 


T> 

MODEL 

68 ,04- 

80.10 

86,97 

107.13 

115.1^ 

^HP 

PLANT 

72.00 

82.00 

86,07 

1C4-.00 

114-.00 

^HP 

MODEL 

^99 #60 

504-.2O 

505.75 

518,70 

521.50 
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The steady state values of the state' variables 
obtained from the model ^-t 90 MW are listed below in Table 


if A. 


TABLE if .if 

STEADY STATE VALUES OF STATE VAKEABLES FOR 90 M¥ 


S.No. 

1 STATE 

VARIABLE I 

STEADY STATE VALUE 

1 ' 

'^CB 

Kg/Sec. 

74.00 

2 

^R1 

Kg/m^ 

7.88 

3 

H 

^^ED1 

Koal/Kg 

771.00 

if 

^IM1 

o„ 

C 

460.00 

5 

%2 

Kg/Sec. 

74.00 

6 

^E2 

Kg/m^ 

6.20 

7 

^E02 

Kcal/Kg 

844.00 

8 

^E02 

Kg/Sec. 

74.00 

9 

^IM2 

°C 

590.00 

10 

%P1 

Kg/Sec. 

72.50 

11 

^MP2 

Kg/Sec. 

70,50 

12 

%P3 

Kg/ Sec. 

68.80 

13 

%P0 

Kg/ Sec. 

65.60 

14 

’■^LPO 

Kg/ Sec^ 

61,88 

15 

^CM 


45.00 

16 

^HM1 

°C 

110.00 


Continued 
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S,No. ' 


17 

18 

19 

20 
21 


STATE VARIABLE 


^HM2 

^HM3 

^HM5 

^HM6 


o. 


STELDT STATE VALUE 

143.00 

220.00 
315.00 
4-50.00 
34-8.00 


DISCUSSION 

The comparison of governing stage exit pressure 
determined by model to the plant data is shown in Table 4-. 3 
The difference in the two values is highest at 5*5^ for 
70 MW load. At 80 the difference is 2.4-;^ on the lower 
side. For 90 MW both model and plant values nearly match. 
This is due to the fact that the model constants were deter- 
mined using the steady state values at 90 MV/. At 100 MW 

the difference in model and plant values is at 2.8^ on 

at 110 MW load. 

higher side and there is a mismatch of only 0»7%j/ The 
error increases as the load value moves down from the 
90 MW level. On the higher side this trend is not reflec- 
ted. This can be due to various assumptions, simplifica- 
tions and guesses made wherever the plant physical 
constants were not available. 
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4.3 TEANSIMT RESPONSES 

4-, 3.1 Plant Responses For a Load Drop of 10 MW 
From the Steady State at 90 MW 

A plant test was conducted on the second 110 MW 
unit of the Panki Thermal Power Station, Kanpur. First, 
when the plant was operating at steady state at 90 MW, 
various control loops were put on manual i.e. the plant was 
put on open^loop. Then, the load setting was reduced hy 
10 MW. All the indicators available in the control room 
were observed by a group of 20 persons with one person 
taking down thf' reading for each indicator. The readings 
were taken manually at a time interval of 2 '- sec. up to 
20 sec. then at a inter-val of 5 sec. up to 5 minutes. 

Various recorder charts that were operative in .the control 
room were also obtained. Out of these, few parameters are 
compared with model tests. It can be pointed out here that 
very little instrumentation was available on the turbine side. 
Further, the instrument calibration was not checked before 
the test andwDiViy persons taii'ng do'^n thJ readings • 
of the indicators were not experienced in this sort of work. 
Hence the results obtained from plant test can only act 
as guide line 

4 . 3,2 Model Response For A Load Drop of 10 MW From 
Steady State of 90 MW 

The model responses for a load drop of 10 MW are dis- 
cussed in the following articles. For obtaining these, 
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the total input flow to turbine, was decreased by a step 
to its new value at 80 MW. 


Ei^h Pressure Turbine 

A drop from 90 MW to ^ MW in load results in a 
step decrease in the input flow to the turbine. The res- 
ponses obtained for the outlet flow from H .P , turbine are 
shown in Fig. 10. The outlet flow from the turbine (Wqj^) 
takes time (equal to its time constant) for adjusting to 
the new value of W,p. This has been brought out in Fig. 10a. 
The turbine outlet pressure undergoes a rapid change in 
pressure as is evident from equation 3*16. The response 
obtained for this pressure from the plant tests is compared - 
with the model response in Fig. 10b. As is evident from 
figure the pressure in the plant stabilises at a higher 
value than the model. This can be due to some error in 
instrumentation or in observing the readings of the indica- 
tors due to the inexperience of the observer. The tempera- 
ture of the outlet steam from the H.P. turbine follows the 
same trend as the pressure (Fig. 10c). This happens because 
the temperature at outlet of the turbine depends on the 
inlet pressure and temperature both of which fall down very 
rapidly with a decrease in load. 
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^.3.2.2 Triflex and Exit Reheater Responses 

The responses obtained for the two reheaters are 
shown in Figure 11. Figure 11a shows the variation of 
Triflex reheater outlet flow. The flow, starts going 

down rapidly and by 6 sec. it becomes more or less stable. 
This is evident from equation 3*31 • Because of the large 
volume Involved in a reheater, it behaves like a storage 
tank. The pressure after the first reheater also 

decreases as shown in Figure 11b. This is explained by the 
following argument. is dependent on the inlet pressure 

to the reheater, the flow and the density. The inlet pre- 
ssure to the reheater goes down fast because of the part 
played by the H.P. turbine dynamics. Further, because of 
the drop in flow, there is storage effect in the 

reheater as a result of which the density goes down. Com- 
bined effect of these results in a gradual decrease of 
pressure. Figure 11c shows the response of the heat transfer 
rate from metal to steam (Qp 3 -])» the flow passing through 
the re heater goes down the heat added to the steam from the 
metal starts going down due to reduction in convective heat 
transfer as a result of decrease in steam velocity, Ibis 
gives rise to two effects. First the heat stored in the 
metal walls goes up which results in an increase in the 
reheater metal temperature, second because of the less 
heat being added to the steam now, the temperature of 
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steam Tj^ 3 -] starts going down. Both these effects tend to 
increase the heat transfer re;.te As a result first 

Qj ^31 starts going down because of the dominant effect of 
flow decrease and then this downward trend is slowed down 
due to the effects that the metal temperature and steam 
temperature exercise on Qg 3 ^ • A;rter the flow has stabilised 
around 6 seconds, its effect on the vanishes and 

because the metal temperature and steam temperature take 
little more time to stabilise there is a marginal increase 
in the heat transfer rate Qp 3 ^ • The steam temperature 
follows the trend shown by Qp 3 ^ (Big. 11 d). The steam 
density response has been shown in Fig, lie. 


The Exit reheater outlet flow behaves the same way 
as the Triflex reheater outlet flow. Tliis Exit reheater 
outlet flow first goes down and then stabilises after 7 secs. 
This is shown in Fig, 11f. The outlet pressure of the Exit 
reheater also shows the same response as that of the Triflex 
reheater. In this case response for comparision to plant 
test was available and both have been drawn in Fig, 11g, 

T’l;: plant test response shows a slower decrease than the 
model response up to 1 sec. and then both show a appre- 
ciable fall in the rate of decrease. Exact matching was 
not obtained due to the same reasons as mentioned earlier. 

The metal to steam heat transfer rate Qp 32 responds in a 
way shown in Fig. 11h. The shaip) rate of decrease, unlike 
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the Triflex reheater case, can be attributed to the fact 
that tube metal mass in the Exit reheater is larger than 
the Triflex reheater and the metal temperature increase 
effect due to heat storage in the metal walls, which tends 
to increase the heat transfer rate, goes down. As a result 
there is no restraint on Qj^32 initially. But this decrease 
of Qj^g2 rise to steam temperature decrease which in 

turn tries to increase the heat transfer rate. The 
enthalpy at the exit of second rehcater follows a very 
similar trend as shown in Fig. 11i. 

4-, 3. 2, 3 Medium Pressure Turbine 

The responses obtained for the mediun: pressure 
•turbine are shown in Fig, 12. The flow at the outlet of 
first M.P. turbine responds in a way shown in Fig, 12a, 

This flow decreases with the decrease of inlet flow to the 
turbine. The time constants of the individual M.P. turbine 
are of the order of 1 sec. which is very small compared to 
reheater time constant and as a result the flow response 
%P1 nearly the same as that of the second reheater. 

The pressure at the exit of the first M.P. turbine is affec- 
ted by the inlet pressure to the turbine, the flow passing 
through it and the density of the steam. The pressure 
first goes down because of the effect of the decrease in inlet 
pressure to the medium pressure turbine, but as flow also 
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goes down and there is only a marginal decrease in the 
steam density, the pressure goes up as can he seen from 
equation 3.4-1. Ihe temperature of steam at the exit of 
first M.P. turbine follows a similar trend due to its 
dependence on pressure and also due to the effect of second 
reheater enthalpy variation (Fig. 12c). 

The second M.P, turbine behaves in the same way 
as the first M.P. turbine. Flow response is shown in 
Fig. 12d, There is no noticeable difference in the response 
from the first one, a fact, which is corroborated by the 
low time constant of the M.P. turbine. The response for 
pressure is shown in Fig, 12e and the temperature response 
is shown in Fig, 12f . The trends are similar to the ones 
obtained for the first M.P. turbine. 

The third M.P. turbine also behaves like the first 
two M.P, turbines and responses are plotted in Fig. 12g, 

12h and 12i. Figure 12g shows the outlet flow response. 

It also shows no noticeable difference in the dynamics 
because of the low time constant and the large time scale 
which had to be used for this. Outlet pressure response 
has been shown in Fig, 12h and the outlet temperature res- 
ponse is shown in Fig. 121, 

The fourth M.P. turbine responds in a similar 
way as the first three ones and the outlet flow, pressure 
and temperature responses are given by figures 12j, 12k, 
and 121 respectively. 
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V.3«2.4 Condenser and Feedwater Heaters 

Figures 13a- to I3g show the responses in metal to 
water heat transfer rates for condenserj low pressure feed 
water heater No, I 5 2, 3? ^ a.nd the high pressure feed water 
heater Nos. 1 and 2, respectively. All show a decreasing 
trend which can he explained by the fact that as the load is 
decreased the quantity of steam extracted also goes down. 

As a result less steam condenses on the tubes resulting in 
lesser heat transfer to the feed water passing through these 
tubes. The metal temperatures of the tubes in the feed water 
heaters and condenser do not change much because of the thin 
tube thickness, 

^, 3.3 Model Responses For 20 Percent Increase in 

Heat Transfer Rate In ReheatersFrom Flue Gases 

In another transient study conducted, the heat 
transfer rate to the tube metal walls from the flue gases 
is increased by 20^^ step. The responses obtained for this 
are shown in Fig, 1^, "^■^guros 1^b and l4-c show the 

responses of metal to steam heat transfer rate, steam enthal- 
py and metal temperature for the Triflex reheater. As the 
heat transfer rate from gas to metal is increased by a 
step the tube metal temperature also increases slowly due 
to the large metal mass involved. This increase in metal 
temperature increases the heat transfer rate to the steam 
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from the metal Qj^g^ • ^is is brought out by the Fig. l4a. 
The metal temperature variation is shown in Fig. l4c. The 
increase in heat content of the steam is evident in the 
response of steam enthalpy (Fig. 1^b) , Similar responses 
are obtained for the Exit reheater and are shown in Fig. l^-d, 
l4e and l4f . The Exit reheater metal temperature increases 
more slowly than the Triflex reheater because the metal 
mass in the Exit reheater is approximately two and half 
times that of the Triflex reheater. This trend of slow rise 
in metal temperature is reflected in the metal to steam heat 
transfer rate and the steam enthalpy. Other parameters 
downstream of the reheater do not show substantial variation 
in their values but a marginal increase in downstream enthal- 
pies and power produced by the turbine is noticed (not shown 
in figure ) , 

4 . 3 .^ Future Work Suggested 

In the present study the model of the turbine for 
a 110 MW unit is developed. Model for the boiler of the unit 
was developed earlier C1T1 - The next obvious step would be 
to combine these two models and to incorporate various con- 
trol schemes in the model. Once this is done many meaning- 
ful tests can be conducted on the plant and the responses 
can be compared with the model . Then the model could be 
used to decide the control strategies and their evaluation. 
Further, the complete model can be used for designing higher 
capacity units by scaling it up. 
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APPENDIX - 1 


PLANT DATA 


The following data has heen taken from the plant 

manual s . 


Steam pressure at the inlet of the throttle valve 
Steam temperature at the inlet of throttle valve 
Steam density at the inlet of throttle valve 
Steam enthalpy at the inlet of throttle valve 
Number of throttle valves 
Number of nozzles in the first nozzle ring 
Number of nozzles in the second nozzle ling 
Number of nozzles in the third nozzle ring 
N-umber of nozzles in the fourth nozzle ring 
Number of stages in the high pressure turbine 
Number of stages in the medium pressure turbine 
Number of stages in the low pressure turbine 


= 130 Kg/cm^ 

= 535 

= 37.8 Kg/m^ 

= 819.8 Kcal/k 


= 4 
= 8 
= 8 
= 7 
= 7 
= 8 


= 12 

= 4 on each 
side 


Number of low pressure feed water heaters 
Number of high pressure feed water heaters 
N-umber of tubes in the Triflex reheater 
Number of tubes in tiie Exit reheater 
Heating surface available for Tiiflex 'rdhieater 
Heating surface available for Exit reheater 
Effective steam volume in the Triflex reheater 


= 5 
= 2 
= 160 
= 224 
= 850 m^ 
= 1714 
= 11*6 m 


I 



Effective steam volume in the Exit reheater = 

Tube metal mass in the Triflex reheater = 

Tube metal mass in the Exit reheater = 

Maximum design metal temperature of Triflex 

reheater = 

Maximimi design metal temperature of Exit reheater= 
Outside diameter of Triflex reheater tubes = 

Outside diameter of Exit reheater tubes = 

Triflex reheater tube thickness = 

Exit reheater tube thickness = 

First extraction steam flow at 90 MW load = 

Second extraction steam flow at 90 MW load ;= 

Third extraction steam flow at 90 MW load = 

Fourth extraction steam flo^^r at 90 MW load = 

Fifth extraction steam flow at 90 MV 7 load = 

Sixth extraction steam flow at 90 MW load = 

Seventh extraction steam flow at 90 MW load = 

Eighth extraction steam flow at 90 MW load = 

Efficiency of high pressure turbine = 

Efficiency of medium pressure turbine = 

Efficiency of low pressure turbine = 
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39.1 

o 

m 

10.33 

. Tons 

30.59 

Tons 

4-6 0 


590 


70 mm 

60.3 

mm 

4-. 5 nmi 

6.0 mm 

1 .680 Kg/ sec . 

1,715 Kg/sec. 

2.616 Kg/sec* 

2.82 

Kg/sec . 

1 .86 

Kg/ sec . 

1.358 

0 

CD 

CO 

4 -. 97 

Kg/sec . 

6.38 

Kg/ sec . 

0.85 


0.85 


0 

CO 

• 

0 



Th© following cons'fcan'cs 3 ,ro dotermined for fhe 
model at the time of initialisations 


^0 

= 0.0813 

Clf 

= 0 . 064-9 


= 0.0208 

^6 

= 0 . 02^1 

^7 

= 0 . 04-65 

^8 

= 0.0567 


= 17.1529 

‘^M ¥1 

= 1.8835 

%MW 2 

= 1.2286 

%iW 3 

= 0.99260 


= O.8O689 

%MW 5 

— 0 . 644-09 

%MW 6 

= 1.2703 

^EMS 1 

= 1 . 4864 - 

^EMS 2 

= 2.7102 

^R 1 

= 30.5 

^R 2 

= 12.1 

^1 

= 0.000554-88 

^2 

= 0.00965 

K3 

= 0.05012 

Kl, 

= 0 . 04-0793 


= 0.090966 
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